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Abstract 

We have studied how the lateral interaction affects the electric conductance, the 
second-order current correlation (shot noise), and the third-order one (skewness) of 
a pair of parallel Al atomic wires. The field operator of wave function is introduced 
to calculate the current-current correlations. The corresponding wave functions are 
self-consistently obtained by iteration according to the Lippmann-Schwinger first- 
principles calculation. The results show that when the distance between two wires is 
sufficiently small, the bonding interaction near the Fermi level will induce a spatial 
derealization of electrons. This leads to an electric conduction value greater than 
the sum of two uncorrelated atomic wire, and a side-band peak is generated around 
the maximum conductance. This correlation can also be observed from the shot 
noise and the skewness. In addition, we find that the three-order Fano factor is 
negatively related to the conductance. 
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Introduction 



Molecular transistors have become an interesting and popular research topic 
due to recent advances in nanotechnologyjl]. From macroscopic aspect, the 
electric conductance of conductors with the same structure is linearly pro- 
portional to the amounts of charge carriers. However, such relationship no 
longer sustain because the interference between atomic bonding and orbitals 
is very important when devices are scaled down to the size of molecules. Lately, 
the studies on electric conductance of molecular wires and monolayers have 
been performed by mechanically controllable break j unctions [2j[3j[4j and self- 
assembly on a solid surface [5] [6]. Nonetheless, few have focused on the cor- 
relation of distance between the atomic/molecular chains. As a consequence, 
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this paper aims to investigate the interactions between wires and the current 
correlations of various moments. 

According to Levitov-Lesovik transport formula [7], all moments of carrier dis- 
turbance in the long time limit may be expanded by transmission probability 
T n (with n the channel index), such as Si oc J2 n T n and 52 oc J2 n T n (1 — T n ). 
The former gives the average current and the later defines the fluctuation of 
current (shot noise) at zero temperature. Recently, the third moment S3, i.e., 
the asymmetry of the distribution function of current fluctuation, is of partic- 
ular interest, and has also been experimentally probed by many groups. [8 J [9 J 
Basically, the disturbance originates from the discontinuity of electron itself. 
In Schottky's vacuum tube experiment, S p = 2e(I), and F 2 = S2/S p = 1, 
which presents no correlation existing between electrons. Since carriers cross 
the junction randomly and independently, this electron transfer can be de- 
scribed by the Poisson statistics. In nanostructures and mesoscopic devices, 
however, there will be correlation between electrons in the transport process 
due to the electronic structure, electromagnetic forces, the Pauli exclusion 
principle, and the Coulomb interactions [TU]. As a consequence, its Fano factor 
may be greater(super-Poisson) or less(sub-Poission) than unity. The shot noise 
can serve as a convenient tool for the probing of electronic structure and of 
the charge amount of a unit carrier, helping understand the nature of trans- 
port. Until now, most studies adopt model-based parameters to simulate the 
transport process. Nevertheless, first-principles calculation offers a possibility 
to accurately link electronic structure to all moments of the current correla- 
tions, which is superior to the model calculations. For instance, Chen and Di 
VentrapJ] found that the shot noise of Si atomic wires of different lengths os- 
cillates as a function of atom number in the wire. Bin Wang and Jian WangjB] 
performed first-principles calculations based on the non-equilibrium Green's 
function to obtain the shot noise properties and found that current increases 
almost linearly with the bias voltage for all Pt atomic structures. All these 
calculated Fano factors well coincide with the recent experimental data. 

In the transport process, charge carriers enter the nanoj unction from a lead. 
Accordingly, it can either arrive at the second lead or reflect back. As a pro- 
totypical example, in this paper this junction is made of a pair of atomic alu- 
minum wires that are sandwiched by two planar jellium electrodes. The main 
reason for choosing Al atoms is because its electronic configuration (3S 2 3P 1 ) 
is simpler than Au atoms so that it is more convenient to observe orbital 
bonding and interference of channels (three valence electrons: sp z , p x and 
p y [14j). Besides, studies on aluminum nanojunction have also been exten- 
sively applied to experiments[13j[15j. In the experiment by Van den Borm 
and van Ruitenbeek[16j, the Fano factor of Al contacts was around 0.3~0.6, 
which proves that the electrons within the junction travel in order, instead of 
random motion as in classical physics. 
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In general, transport problems deal with non-equilibrium states, and the count- 
ing statistics of current is often related to the bias voltage between termi- 
nals. Therefore, these issues cannot be tackled directly by the equilibrium 
ground-state of non-interaction electrodes. To overcome this difficulty, we 
adopt the density functional theory (DFT) coupled the Lippmann-Schwinger 
equation (LS-DFT) and perform self-consistent calculation for the ground- 
state wave function of systems with nanojunctions at limited bias voltages. 
This approach is essentially equivalent to the DFT combined with the non- 
equilibrium Green's function [T7](NEGF-DFT). One different is that the wave 
function used in LS-DFT consists of the leads and the junction, while the 
one for NEGF-DFT is restricted in the junction and a self-energy has to be 
considered in the Green's function[18j. When calculating the shot noise and 
skewness through the field operator of wave function, we are faced with the 
products of field operators and the Bloch-de-Dominicis theory is applied to 
recast these multiple operators into pairs. The results resemble the previous 
study on multi-current correlation functions by the scattering theory. [19] 

In mesoscopic systems, the size of junction is typically greater than the de- 
phasing length of electrons. The wave functions of electrons will lose coherence 
after multiple scattering. The interference of transverse quantum channels is 
neglectably small. However, for the junction of atomic size, the transport of 
electrons can be coherent. Hence the overlap integral of wave function at differ- 
ent transverse momenta K must be kept in describing the current correlation 
of different moment based on wave functions. In other words, the effects of 
multi-channel mixing is included in our calculations. Until now, the analysis of 
current correlations of atom-sized parallel Al chains by first-principles hasn't 
been proposed. To compele this goal, we have developed a current operator 
composed of DFT self-consistent wave functions. Even though the current of 
higher moments is hard to measure and calculate, it can be demonstrated to 
provide a way to obtain the detailed interference effect caused by wave func- 
tions at atomic scale. In the following, we expand the current correlation of 
different moments in the basis of wave functions. The space parameters of the 
two metal contacts connecting parallel Al wires are 5.8 a.u. between Al atoms, 
the Al radius is r s = 2. a.u, and the last Al atom is placed 2a. u. outside the 
boundary of jellium model positive background. 



1 Theoretical Model and Calculation Method 

The Hamiltonian for the whole system is H = Hq + V(r), where Hq = 
J2 -Ek ,,a^K h^Ikh stands for the bare electrodes described by the jellium 

B,K||,a 11 11 

model and V (r) is the scattering potential of the central nanostructure. The 
system's ground-state wave function can be obtained through self-consistent 
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calculation with the Lippman-Schwinger equation, which is. 

*S§ (r) = (r) + / d 3 r, J dh 2 G% (r, r,) V (r u r 2 ) (r 2 ) (1) 



Here, ( r ) represents the unperturbed wave function of the electrodes 

before scattering the nanostructured object, and it can be expressed 

(r) = (27r)~ 1 e* K H' R • ^kk,! where ttg5 (z) is the wave propagating 
along the z-direction. When z — > ±00, the boundary conditions satisfied by 
wave functions on the left and the right sides are 
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T L e ikLZ ,z > 00 

+ R R e ikLZ ,z> 00 
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T R e~ ikLZ 



(2) 
(3) 



The subscript Ky is the momentum of electron with its direction parallel to the 
surface of electrodes, while e tkaZ is the free propagating wave perpendicular to 

it. The electron's momentum is given by k a = ^2m (E — K|| 2 /2m) — v e ff (±00), 

v e ff (z) denotes the energy of the bottom conduction band in electrodes. It 
should be noted that G° E is the Green's function without the addition of the 
nanostructured object. The scattering potential is given by 



v (n, r 2 ) = v ps (n, r 2 )+|(K c [p (n)] - v; c [ Po (n)]) + f dh 3 ^^5 (n - r 2 ) . 

(4) 

V^ s (ri,r 2 ) is a non-local pseudo-potential formed by covalent electrons with 
the core electron being transferred. V xc is the exchange-correlation potential 
derived by LDA, and the final potential is the Hatree potential. p (r) and 
p (r) represent the charge density for the electrodes and the whole system, re- 
spectively. Sp (r) is their difference. In numerical calculation, the box includes 
the middle molecular system and partial electrodes. It is chosen sufficiently 
large so that outer electrodes are not affected by molecules. Within the box, 
2880 plane waves are utilized to calculate this coherent quantum system. 



1.1 First-principles current calculation 
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At first we define a field operator, ^ (r,i) = a EK H (*) *£Kn ( r )' to de- 

B,K||,q 11 11 

scribe the many-body state of the system, is an operator meaning the an- 
nihilation/creation of electron at time t and at electrode a, and it satisfies the 
anti-commutation relations, ja^K^ a £K 2 } = <W<^ (E\ — E 2 ) 5 (K x — K 2 ) f%. 

~ ^l(r)) I^bT\ } is the probability of occupying the a 
electrode, which satisfies the Fermi-Dirac distribution. Based on the definition 
of electric current, we derive that 

i(z,t) = ^J dr ± J dE 1 J dE 2 J dKt J dKae'^-^a^a^i^^ 

(5) 

where /|f Kl , £2 K 2 (r) = (^kJV*^ ~ V(*I iKi )**§ 2 k 2 - At zero tern- 
perature, the average current is given by 

( f (*)) = ^ l dE S ^ l dKi EK,EK (T ± ) , (6) 

1 1 Lb J <J J 

and it is noticed that at steady state, (l(z)) is independent of z. 



1.2 Second moment of the current 



We can extend the derivation and obtain the second and the third moment of 
current. The power spectrum of shot noise has the form 

S 2 (u, T, Zl , z 2 ) = 2nh [d (t, - t 2 ) e 1 ^-^ (SI (z u t x ) 51 (z 2 , t 2 )) , (7) 



with 51 (t) = /(£) — (^1) . When calculating the shot noise, we have to deal with 
the quantum statistical expectation value of the products (o^a^a^a^). 



With the help of Wick-Bloch-de Dominicis theorem, it can be decomposed as 

( a E^K 1 a kK 2 a Ek i K3 a E4-K4) = ( a l^K 1 a B 2 K 2 ) ( a E^K 1 a E2K 2 ) + ( a EtK 1 a E4K4) ( a £ 2 K 2 a i?tK3 

Substituting this result into Eq. (JJJ), we have 



S 2 («, T; Zl , z 2 ) = 2nh (—) £ / dEf E+nui (l - /|) / dr 1± [ dr 2 {8) 

l / ^ K - 2 ^+^K 1 ,BK 2 ( r l) IeKi,E+KuK 2 ( r 2) • 
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Because shot noise reflects the current correlation at zero temperature, the 
steady-state (cj — > 0) shot noise is thus found as 



(eh \ Efr 
■) J2 / rfr l-L / dr 2± 
m V a ,/3=L.R JE ^ J J 

dKi I dK 2 lE+nwK 1 ,EK2 ( r l) I EK 1 ,E+hLoK 2 ( r 2) 



(9) 



1.3 Third moment of the currenM 



With similar approach, the third moment of the current is defined as 



5*3 (uj,uj',z 1 ,z 2 ,z 3 ) = 2Tth(^^ J d(ti-t 3 ) J d(t 2 -t. 



(10) 



Substituting current difference 51 into Eq. f lTUj) and considering the steady- 
state current (u = u' = 0) at zero temperature (T = 0), we have 



2 ( eh \ r E FR 



S 3 ( Zl ,z 2} z 3 ) = (2nhy^—j J e dE j dr 1± j dr, ± j <lr, ± j r/K, / dK, j rf/H : ) 



x 



}RL 



(r 2 ) I 



LR 



'■EK,.IK 2 ( r l) IeK 2 ,EK:; 'EKa.EK 



x (ra) - 



^EK!,EK 2 ( r l) ^K!,£;K2,£;K3 ( r 2) ^BK 3 ,EKi ( r 2 

Note that the current density in S 2 and S 3 keeps interference between wave 
functions with different transverse momenta K, which does not coincide with 
Eq. (©. 



2 Results and Discussion 



To analyze the conductance, the second-moment shot noise, and the third- 
moment skewness, we have derived formulas for different moments of the cur- 
rent correlation through the first-principles DFT combined with field oper- 
ator of wave function. That is, Eq.([6]), Eq.([9]), and Eq. (fTT|) . The Lippmann- 
Schwinger equation is utilized to self-consistently calculate the scattering wave 
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Fig. 1. The differential conductance of two parallel four-atom wires as a function 
of wire separation. The dashed line stands for the sum of conductance of the two 
isolating atom chains. 

function of the whole system. Figure [1] shows the conductance in the unit of 
Gq[= e 2 /h] versus the distance between two parallel wires. 

We find that the differential conductance has a maximum value of 3.68 GO 
at the distance of 5 a.u. and a side peak at the distance of 5.5 a.u. As the 
distance is increased, the conductance still presents oscillating changes (after 
8 a.u). At last, the conventional conductance is calculated when the distance 
between two wires is effectively infinite, i.e., two times the conductance of a 
single wire. In order to investigate the modification on the conductance by 
the interactions between wires, we study the distribution of carrier density at 
d = 4.5~7.5 a.u. 



From Fig. [2j it can be seen that the cr-bonding state is formed through the 
interactions between the S and the P z orbitals in the atomic wire comprising 
Al atoms (3S 2 3P 1 ), which decides the distribution of charge carrier between 
wires. While the distribution along Al wires is determined by the IT-like in- 
teraction between the P x and the P y orbitals. A comparison of carrier density 
distribution shows that the minimum density of 0.23 occurs at d=5.0 a.u. and 
that the minimum of 0.16 also takes place along the wire. It indicates that 
the derealization of carriers reaches its maximum at this distance. The de- 
localized electrons do not belong to any bonds or atoms, and few electrons 
accumulate in this region. Finally, large amounts of electrons rapidly travel 
through the system and arrive at the other electrode, causing the maximum 
current. We observe the distribution of DOS for further understanding the re- 
lationship between derealization and distance as in Fig. |3l Note that a small 
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Fig. 2. The spatial distribution of carrier density 

bias voltage is exerted, , and the Fermi energy in the left electrode is set to 
be zero. 

Figure [3] shows the DOS(Ep) of parallel atomic chains at varied distances. 
To observe its local density of states in the atomic wires, the projected DOS 
is utilized. Take the case of d=4.5 a.u. for example, we understand that the 
first two peaks (located between E = — 0.15~ — 4.5) are mainly contributed 
by the S and P z orbitals. Nevertheless, the orbital of P xy influences the region 
near the Fermi energy. When increasing the distance, the bump between the 
first two peaks disappears gradually and eventually become one single peak. 
It is noticed that the peak shape of P xy orbital is modified and the extent has 
positive relationship with the change in the DOS with the distance (top-left 
inset of Fig. [3]), but the peak position is not shifted. It means that the change 
in carrier density caused by the distance for cx-type interaction is far greater 
than that for Il-type interaction. 

Normally, the measured conductance corresponds to the change of the elec- 
tronic structure at the Fermi energy. In this case, it reflects the exact change 
of P xy orbital. Comparing the top-left inset of Fig. |3]to Fig. [TJ we find that the 
DOS (Ep) cannot completely decide the tendency of conductance. For exam- 
ple, the DOS(-Ep) exhibits a parabolic curve at 5.5 to 7 a.u., forms a plateau 
at 7 to 7.5 a.u., and then declines linearly. Such results disagree with those 
for conductance, which shows linear decrease only. This owes to the fact that 
electron transport is a dynamic process. In calculating conductance, we need 
to consider the transition rate of electrons in addition to the DOS. The tran- 
sition rate versus the distance in the nanojunction is shown in the top-right 
inset of Fig. [3j It is observed that the electronic transition rate has its max- 
imum at d — 5.2 a.u., shows a valley at 7.5 a.u., and climbs up afterwards. 
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Fig. 3. The distribution of DOS in the AL wire between two electrodes at d=4.5, 5, 
and 5.5 a.u. The top insets present the DOS of single wires, i.e., the distance is 
sufficiently large. The Fermi energy on the left side is zero, and the bias voltage is 
O.OleV. The gray line denotes the transport regime. The top-left inset gives the DOS 
at the Fermi level as a function of wire distance. The dashed bar on the right shows 
the value for infinite distance. The top-left inset is the tunneling rate of electrons 
near the Fermi level. 

So, when the DOS (E F ) is multiplied by the tunneling rate, a side-band peak 
exists near the maximum derealization distance, resulting in a similar curve 
as conductance. 

To realize the relationship between moments of counting statistics of current 
fluctuations in the Aluminum chains, we focus on how the second-order Fano 
factor (related to the shot noise) and the third-order Fano factor (related to 
the asymmetry of the distribution) change with the distance. In Fig. HJ F 2 first 
decreases and then increases, and two dips are respectively formed at d — 5 
and 5.5 contrary behavior with respect to the conductance (two peaks). 

According to Ref. [21], the suppression of shot noise in metallic conductors 
arises from open quantum channels. Here we know that the channels are 
mainly from the Pxy orbital, especially at d — 5 and 5.5 a.u., and that the 
interaction between p orbitals further induces a spatial delocalization of elec- 
trons, leading to the sudden rise of current. The bottom inset of Fig. 4 plots 
the third moment versus the separation distance. It is found that F 3 is weakly 
negative at varied distances, indicating the distribution function of the num- 
ber of electrons transferred through a conductor within a given time period is 
negatively skewed, and it is also strongly negatively correlated with the con- 
ductance. Different from the studies on single C4 or C5 chain, where F 3 and 
conductance are always positively correlated. In this case, we find that F 3 is 
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Fig. 4. The second- and third-order Fano factors versus the distance between two 
parallel Al chains. The right dashed line denotes the value at infinite distance and 
at bias voltage of O.OleV. 

positively correlated before d =7.5 a.u. and reveals highly negative correlation 
at greater distances. It should be noticed that the least conductance occurs 
at d =8 a.u. while the F3 value is strongly skewed, implying it is difficult for 
charge carriers to pass through the junction. 
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